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Abstract 



A 9 term, which couples to topological charge, is added to the two- 
dimensional lattice CP 3 model and U(l) gauge theory. Monte Carlo 
simulations are performed and compared to strong-coupling character 
expansions. In certain instances, a flattening behavior occurs in the 
free-energy at sufficiently large 9, but the effect is an artifact of the 
simulation methods. 
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I. Introduction 

Formally, the inclusion of a 9 term in a theory does not affect the equations 
of motion. For this reason, it was once thought that 9 terms were irrelevant. In 
fact, this appears to be the case for the four-dimension U(l) gauge theory. Adding 
the 9 term / d^xF^F^ix) to the action is believed to not affect physics. The sit- 
uation for a non-abelian theory is different. The importance of the 9 term Sg = 
g 2 9 J d 4: xF^ u F^ u (x)/(327i 2 ) was realized when instantons were uncovered in four- 
dimensional Yang-Mills theories Instantons represent barrier-penetration pro- 
cesses between different classical n- vacua. An n vacuum \n) is obtained from the 
perturbative vacuum through a non-trivial gauge transformation carrying winding 
number n. Instanton-barrier-penetration effects imply that the true vacua are linear 
combinations of n vacua. These vacua \9) are called 9- vacua || |3| and are given 
by \9 >= X^-oo ex P (in9) \n). In a functional integral, #-vacua are incorporated by 
adding the term Sg to the action. Since Sg breaks parity, time-reversal invariance and 
CP symmetry when 9 ^ or 9 ^ ir, the strong interactions explicitly violate these 
symmetries for < 9 < n. In fact, not only do instantons contribute to Sg but all 
topological quantum fluctuations also do so. Such topological fluctuations are known 
to exist because they contribute significantly to the rj' mass 

In QCD, phases in the quark mass matrix M. also contribute to CP violation. 
However, only the combination 9 e ff = 9 + ArgDetAi is relevant for strong CP viola- 
tion for the following reason. Through redefinitions of quark fields, one can eliminate 
the CP violating phases in the quark mass matrix. However, to eliminate one CP- 
violating phase, a (7(1) axial rotation is used. Via the axial anomaly, ArgDetM. 
re-emerges as a coefficient of Sg. Hence, the physical effective theta parameter is 
9 e ff. The strongest constraint on 9 e ff comes from the electric dipole moment of the 
neutron. Compatibility with experimental bounds requires 9 e ff < 10~ 9 [0, The 
strong CP problem in QCD can then be phrased as the question of how 9 e ff can 
naturally be so small. The strong CP problem is thus a fine-tuning issue. 

In a pure Yang-Mills theory, the strong CP problem involves vacuum dynamics. 
Vacuum physics is related to the long-distance behavior of a theory and hence non- 
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perturbative and strong-coupling effects. For four-dimensional Yang-Mills theories, 
this regime is not well understood. Thus it is useful to consider simpler systems, 
such as the two-dimensional CP N ~ l models. They have features in common with 
four-dimensional Yang-Mills theories: They are asymptotically free, possess instanton 
solutions, and have 9 vacua. 

Insight into the CP N ^ models comes from strong coupling expansions JT5 



Monte Carlo simulations and the use of the large N limit [^| [HJ. As N — ► 

oo, the CP N ~ X model becomes a system of N free particles and N free antiparticles, 
and there is no 9 dependence: #-vacua are degenerate in energy. However, the first 
1/N correction lifts the degeneracy of the 9- vacua, with the energy separation between 
9- vacua being of order 1/N. The first 1/N correction also leads to a quantum- 
mechanically generated £7(1) gauge field. This gauge field produces a linear potential 
between particles and antiparticles. Hence the system confines, although the string 
tension, in this approximation, is again of order 1/N. Actually, it is expected that 
the confining force becomes stronger and stronger as higher-order 1/N corrections 
are included since the CP N ~ l models exhibit superconfinement ||17|| . The dramatic 
changes in behavior in going from leading order to next-to-leading order are indicative 
of the singular nature of the 1/N expansion in the CP N ~ l models. The singular nature 
is also reflected in the fact that strong coupling and 1/N expansions do not commute 



Tq , [17], For this reason, results from Monte Carlo methods should probably be 



trusted over results from 1/N methods. 

In the presence of a 9 term, G. Schierholz and co-workers in pioneering work have 
performed Monte Carlo simulations on the CP 3 model. Results for the free energy 
were interesting and unexpected. For fixed inverse coupling (3, a dramatic change in 
the free energy behavior occurred at a critical value 9 C of 9. The free energy / per 
unit volume was well represented by 

j a((3) 9 2 9 <9 C 
1 \c(P) 9>9 C ■ [L - L) 

Hence, for 9 > 9 C , the free energy had no dependence on 9. In two dimensions, the 
string tension cr(e,9) for external particles of charge e in a theta vacuum \9) can be 
computed from the free energy using <r(e, 9) = f(9+2ire) — /(0).p2fl Hence, for 9 > 9 C , 
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the string tension vanishes for particles of sufficiently small charge. In other words, 



confinement for small external charges is lost. Furthermore, the simulations in p8 



indicated that 9 C went to zero in the continuum limit in which the coupling g goes to 
zero. 



In strong coupling character expansions of the free energy were obtained. 
Near infinite coupling, no flattening behavior of the free energy like that of Eq. ( |1 . 1| ) 
was seen. However, at smaller g, a peak in the free energy occurred. Although 
higher corrections might change this behavior, the strong coupling series provided 
some support for the form of the free energy given in Eq. ( |1 . 1|) . 

Many two-dimensional systems with a U(l) gauge field, such as the Schwinger 
model, possess a cusp in the free energy at 9 = it. The cusp signals the spontaneous 
breaking of CP invariance.||33|| At infinite coupling, the lattice CP N ~ l models also 
undergo spontaneous CP breaking at 9 = 7T.||18|| Hence, a phase transition in the 
(jpN-i moc i e } s at C = 7T is not unexpected. What is interesting about the numerical 
studies [28] of the CP 3 model is that, for sufficiently weak coupling, 9 C moves away 
from Ti and decreases with decreasing g. This picture suggests that, to obtain a 
continuum confining theory from the lattice CP 3 model, one must tune 9 to zero. 
This then also suggests that continuum confining CP N ~ l models must have 9 = 0.Q 
If the analog of this statement were true for a four- dimensional Yang-Mills theory, 



then the strong CP problem would be solved. In fact, preliminary studies J34| of 
the free energy of the four-dimensional SU(2) Yang-Mills theory show a free energy 
behavior similar to the one in equation (|1 . 1|) . However, in four-dimensional Yang- 
Mills theories, the string tension is not related to the free energy /. Hence, the 
analog of the argument for the CP 3 model, namely, that 9 must be tuned to zero to 
obtain confinement, does not necessarily hold. Nonetheless, if 9 must be less than 
9 C for some other physical reason, and if 9 C goes to zero as g — > 0, then the strong 
CP problem would be solved in the Yang-Mills theories. This is interesting because 
it has been suggested that the strong CP problem might be solved naturally within 



the pure Yang-Mills theory. [35] In fact, such a result occurs in a 2 + 1 dimensional 
model: the Yang-Mills sector generates a relaxation field, which acts like the axion 



-'^One would still need to explain why confinement is necessary for the continuum limit. 
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in the Peccei-Quinn mechanism p6l 137], |38|. In ref.[p5|, criteria were established to 



determine when such a natural relaxation mechanism arises for theories in arbitrary 
dimensions. It has not been yet possible to determine whether these criteria are 
satisfied for four- dimensional Yang-Mills theories. 



One unusual feature of the Monte Carlo simulations of ref. [28| is the volume 
dependence on the free energy for sufficiently large 9. As the volume V increases, 9 C 
decreases. This leads to a situation in which the free energy differs by large factors as 
V varies. For example, at (3 = 2.7 and 9 near 7r, the free energy on a 46 x 46 lattice 
was three times the free energy on a 72 x 72 lattice. See Figure 3 of ref. ||28|| . Since 



quantities, like the free energy, should quickly converge in the thermodynamic limit, 
there are only three possible explanations: 

(i) For sufficiently large 9, there are light, perhaps massless, modes in the system 
which cause finite-size effects. 

(ii) A systematic effect occurs which leads to numerical results that differ from true 
results. 

(iii) There is some unknown explanation not covered by (i) or (ii). 

One purpose of the current work is to try to determine whether (i) or (ii) occurs. 
We perform Monte Carlo simulations on an exactly solvable U(l) lattice gauge theory. 
By comparing numerical and analytic results, much insight into simulating systems 
with a 9 term is gained. Our results are presented in Section III. We will argue that 
systematic effects can lead to anomalous flattening behavior in the free energy, as 
described by Eq. ( |1.1| ) for the U(l) lattice gauge theory. In Section II, a general 
analysis of simulating systems with a 9 term is presented. Section II provides a 
mechanism by which anomalous flattening behavior can arise. In Section IV, a lattice 
CP 3 model is treated. By comparing Monte Carlo data with analytic strong coupling 



series from ref. ||19|| , flattening behavior can be shown to be anomalous for at least 
two simulations. Some additional results and remarks are presented in Section V. In 
Section VI, a summary and final discussion is given. 
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II. General Issues Concerning Simulations with 6 Terms 

In this section, we discuss some issues concerning simulations of an arbitrary 
lattice theory in the presence of a 9 term. In particular, a general error analysis can 
be carried out. From this analysis, one concludes that, for sufficiently large volumes, 
a limiting 9 exists, beyond which reliable measurements of the free energy cannot be 
made. 

For a fixed V, let P(Q) be the probability of having a configuration with topologi- 
cal charge Q in some system. Let Pmc(Q) be the corresponding quantity as measured 
in a Monte Carlo simulation. Assume that the Monte Carlo updating procedure gen- 
erates configurations proportional to Boltzmann weights. Below, this assumption is 
relaxed. If N MC {Q) is the number of times that configurations with topological charge 
Q arise in such a simulation, then 

In a typical simulation, the measured Pmc{Q) differ from the exact P(Q) by small 
errors 8P(Q): 

Pmc(Q)=P(Q) + SP(Q) . (2.2) 

With enough measurements, \5P(Q)\ <C 1. In most systems and the ones considered 
in this work, P(Q) falls off with Q, so that P(Q) > P{Q') for \Q\ < \Q'\. A criterion 
for a simulation to have good statistics is that 5P(Q) <C -P(O). 

Let f(9) be the difference between the free energy T{9) of a system with a 9 term 
and the free energy of a system with 9 = 0: 

f(9)=F(9)-F(0) . (2.3) 

Typically, f(9) is an increasing function of 9 for < 9 < it. The free energy difference 
f{9) is constructed from P(Q) using 

exp(-V7(0)) = $>(Q)exp(i0Q) • (2.4) 

Q 

Normally P(-Q) = P(Q), so that f{-9) = f(0). 
2 Evcntually the falloff is rapid. 
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In a Monte Carlo simulation, an approximation /mc($) to f(9) is obtained by 
using Pmc(Q) in lieu of P{Q): 

exp (-Vf MC (9)) = Y j Pmc{Q) exp (i'0Q) = exp (-Vf(9)) + , (2.5) 
Q 

where 

5Z(0) = 6P (Q) ex P W) • ( 2 - 6 ) 
Q 

Hence, 

- V/ wc (0)) = log [exp {-Vf{9)) + 5Z{9)\ . (2.7) 

Since /(#) is an increasing function of 9, an accurate measurement of f(9) for < 
9 < 9q is obtained if 

\8Z(9)\<exp(-Vf(9 )) . (2.8) 

In particular, since /(0) = and <C 1, there is always a region near # = 

for which f{9) can be measured in a Monte Carlo simulation. However, away from 
9 = 0, Eq. (|2.8|) implies that accurate results are obtained only if the error 5Z(8) 
is exponentially small in V. If Eq. ( |2.8| ) is satisfied with 9q = it, then a reliable 
measurement of f(9) can be made throughout the entire fundamental region < 9 < 

7T. 

If the inequality in Eq. ( |2.8| ) is not satisfied for some 8q, then one of several 
possibilities may arise. If, on one hand, exp (— V f(9 ))+5Z(9 ) < then the argument 
in the log on the right-hand-side of Eq. Q2.7| ) becomes negative and one will not be able 
to extract /(#o) from the measurements of the probabilities Pmc{Q)- 111 Monte Carlo 
simulations for such a situation, a growth in the errors of fMc{9) will be observed 
as 9 approaches 9 , and /mc(^) will eventually not be measurable. If, on the other 
hand, exp (— Vf(9 )) + SZ(9 ) > then fMc{9) will be measurable at 9 = 9 but the 
results will not be accurate. 

For sufficiently large V, a value of 9 exists beyond which it is impossible to reliable 
compute f{9). If we call this value 9b, then 9j, is the maximum value of #o for which Eq. 
( P-8|) is satisfied. The value of 9b depends on the statistical accuracy of the simulation. 
As V gets larger, 9b decreases unless an enormously large number of measurements are 
undertaken to reduce statistical errors. For large V, obtaining enough measurements 
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becomes, in any practical sense, impossible. Clearly, it is more difficult to measure 
f{6) throughout the entire fundamental region of 9, as V gets larger. 

It turns out that in most Monte Carlo simulations, there is a tendency for 

\6P(0)\ > \6P{1)\ > \SP{2)\ > ... . (2.9) 

The reason for Eq. ( |2.9| ) is explained in the next paragraph. For the sake of argument, 
suppose that |5P(0)| is much larger than \5P(Q)\ for \Q\ > 1. Then, from Eq. (|2.8|) , 
one concludes that 

/(ft)«~|log|M>(0)|| . (2.10) 
Since Monte Carlo results are reliable for 9 < 9 b , 

fMc(9)^f(9) foi9<9 b . (2.11) 

If, in addition, <5P(0) > 0, then one will findQ 

f MC (9)--^\og6P(0) for9>9 b , (2.12) 

so that a constant "flat" behavior in fMc{9) will be observed. Although one might 
expect the statistical error in fMc{9) to be the order of fMc{9) for 9 > 9 b , we have 
observed in some simulations that the error, as computed by the jackknife method, 
remains anomalously small. This happens in both the two-dimensional U(l) gauge 
theory (see Section III) and the CP 3 model (see Section IV). When 5P(0) > 0, 
fMc(9) will smoothly interpolate between the behaviors in Eqs. ( |2.11|) and (|2.12| ) in 



the region around 9 ~ 9 b . The net result is a form for the free energy which resembles 
Eq. ( |1 . 1| ) . If an error SP(Q) for \Q\ > is almost as large as <5P(0), then the constant 
behavior in Eq. ( |2.12| ) is replaced by a "slightly wavy" almost flat curve. The above 
discussion applied to the case <5P(0) > 0. If, on the other hand, <5P(0) < 0, then it 
will be impossible to measure f(9) for 9 > 9 b . 

There are various ways to measure Pmc{Q)- The naive method is to generate 
configurations proportional to their Boltzmann factor and simply count the number 



Note that since SP(0) <C 1, the right-hand side of Eq. (2.12) is positive 
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Nmc(Q) of configurations with topological charge Q, and then use Eq. (|2.1|). Since 



the statistical uncertainty SN MC (Q) in N MC (Q) is proportional to jN M c(Q) 



\5Pmc(Q)\~c\JPmc(Q) , (2.13) 

for some constant c|] Because P(Q) is a monotonically decreasing function of Q for 
Q > 0, the inequalities in Eq. ( |2.9|) tend to be satisfied. In any particular simulation, 
however, statistical fluctuations may violate Eq. (|2.9|) , particularly if P(Q) is only 
slightly larger than P(Q + 

One commonly used Monte Carlo technique p3| , f£T| is to perform simulations in 
Q-bins. The topological charges Q are restricted to be in a bin of a certain bin- 
size. The bins must overlap - usually one takes the bins to overlap for only one Q 
value. For example, the bins can be taken to be < Q < AQ, AQ < Q < 2AQ, 
2AQ < Q < 3AQ, etc., where AQ is a small positive integer. Here, the bin-size is 
AQ+ 1. The relative probabilities within a Q-bin are measured by generating config- 
urations which are proportional to their Boltzmann factors but with their topological 
charges contained in a Q-bin. Overall probabilities are determined from the relative 
probabilities by matching results at the overlap endpoints of two Q-bins. Compared 
to the naive method, one generates more configurations with large Q configurations. 
Hence, errors in SP(Q) are relatively reduced for large Q. This tends to enhance the 
inequalities in Eq. (fOl). 

A benefit of the binning technique is that it enables P(Q) to be accurately mea- 
sured for large Q, even where P{Q) is extremely small. The effect can be spectacular: 
often P(Q) which are orders and orders of magnitude smaller that P{0) can be mea- 
sured. 

Another simulation method introduces a weight factor w(Q). One then generates 
Monte Carlo configurations which are proportional to (the Boltzmann factor) x 
w(Q). If Nmc{Q) is the number of configurations with topological charge Q which 
are generated by such a procedure, then Eq. (|2.1| ) is replaced by 

n mc (Q) MQT 1 

Eq' N mc (Q') [w(Q')Y 



4 There might be a weak dependence of c on Q. 

5 This tends to happen when the volume is large and for Q small. 
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Although any weight w(Q) may be used, the trial-probability-distribution method 
tries to choose w(Q) so that all Q sectors are "visited" roughly the same number of 
times. In other words, one guesses a trial probability Po(Q), which might accu- 
rately approximate the true probability distribution P(Q). One then uses w(Q) = 
1/Pq{Q). If one were to pick P so that Po(Q) = P(Q), then a constant distribu- 
tion in Q for Nmc(Q) would be generated (up to statistical fluctuations) and Eq. 
( ^.14j ) would lead to Pmc(Q) — Pq(Q) — P(Q)- Like the binning method, the trial- 
probability-distribution method generates more configurations in the large Q sectors, 
thereby enabling one to better measure P(Q) for Q large. This method also enhances 
the inequalities in Eq. (|2.9|) . 

Finally, one can combine the above two methods by using a trial probability dis- 
tribution Pq{Q) with a Q-binning. This trial-probability-binning method is excellent 
for measuring P(Q) for large Q. However, of the four methods mentioned here, it 
enhances the inequalities in Eq. fl2.9|) the most. 

In summary, current methods of simulating systems with 6 terms tend to generate 
errors for P(Q) which are ordered as in Eq. ( |2.9| ). 

None of the methods solve the 0& barrier problem discussed above between Eqs. 
and ( pj.l^j ): Assuming that 5P(0) dominates in the error 5Z(8), one sees that 



increasing the statistics hardly changes 6b because of the logarithm and the volume 
factor in Eq. ( [2. 10 ). The same statement holds even when SP(0) is not the dominant 



error. One concludes that all simulation methods must fail for sufficiently large 6, if 
V of the system is large. Cluster algorithms might help a little in this regard. [43| 



From the above discussions, the following general guidelines concerning Monte 
Carlo simulations of systems with a 6 term are obtained: 

(1) When the volume is sufficiently big, a limiting 6b arises. For 6 > 6b, the free 
energy cannot be reliably measured. 

(2) As long as finite-size effects are under control, that is £ < V^^ d \ small- volume 
results for the measurement of f(6) are more reliable than large- volume results. 
Here, £ is the correlation length and d is the number of dimensions of the system. 



9 



(3) If a flattening behavior of the free energy f(9) for large 9 is observed, one 
should be cautious that the result is spurious. In particular, one should try to 
see whether |5P(0)| is bigger than the other \SP(Q)\. 



(4) When P M c(Q) is less than \SP(0)\, the contribution to P M c(Q) in Eq. (gT 
need not be included. The reason for this statement is that the contribution of 
Pmc(Q) is lost in the "noise" of the error term 5Z in Eq. (|2.5|) . 



(5) Although the Of, barrier cannot be overcome, Monte Carlo procedures should 
emphasize accurate measurements of P(Q) for Q near 0. In particular, the 
optimal procedure is one for which all \5P(Q)\ are approximately equal. This 
minimizes the chances for anomalous flat behavior in f(9). 

Combining (2) and (3), one obtains another guideline: 

(6) If a large- volume simulation shows a flattening effect for f(6) for 9 sufficiently 
large, but a smaller-volume simulation does not, one should probably trust the 
smaller- volume result. 

Surprisingly, many previous Monte Carlo studies of systems with 9 terms measure 
P(Q) over many orders of magnitude. Point (4) above says that, although this is 
not harmful, it is an inefficient use of computer time if one is interested in measuring 
the free energy. Likewise, point (5) implies that the standard binning and trial- 
probability methods are not optimal because they enhance the inequalities in Eq. 
( [2.9| ). The binning method can be improved by doing more measurements in the Q- 
bin containing Q = and less measurements in large-Q Q-bins. Another improvement 
is as follows. One can adjust w(Q) in Eq. ( |2.14j ) so that the errors \SP(Q)\ are equal 



(up to statistical effects). It is not hard to show that the optimal w{Q) is P{Q). 
This corresponds to using a trial-probability function which is the inverse of the true 
probability. In other words, the optimal weighting method is the antithesis of the 
trial-probability method and the antithesis of what is commonly employed in Monte 
Carlo experiments. 

III. The Lattice U(l) Gauge Theory 



10 



The two-dimensional lattice U(l) gauge theory serves an important test case be- 
cause computer simulations can be compared to the exact analytic results below. For 
related Monte Carlo investigations, see refs. [[IT], (PJ. 

In the lattice formulation of a gauge theory, one assigns an element of the group 
to each link of the lattice. For the U(l) case, such an element is a phase. The lattice 
U(l) gauge theory action is 

s u ^ = (3J2{u P + u;) , (3.i) 

v 

where U p is the product of the U(l) link phases around the plaquette p and where (3 
is the inverse coupling. 

Define a local topological density v v via v p = log (U p ) / (2%), where — tt < log (U p ) < 
ii. The total topological charge Q is given by Q = J2 p v p . The theta term Sq term is 
i9Q, that is, [|5j 

^ tC rm = ^E 1 °g( f/ P ) " ( 3 ' 2 ) 
271 p 

Eq. ( p.2|) is the lattice analog of the continuum #-term action i^- J d 2 xF. 

In d = 2 dimensions, gauge theories are exactly solvable even in the presence of a 
9 term. For periodic boundary conditions, in which case Q is quantized to an integer, 



the result is ||T 



Z(6,(3,V)= £ [z(6 + 2nm,P))} v , (3.3) 

m=— oo 

where 



& f exp(%g)exp(2/3cos( / p )) 
exp(2/3cos(/ p )) 

Here, f p can be thought of as the field strength for a single plaquette: U p = exp(if p ). 
In the infinite volume limit, the free-energy difference per unit volume f(8) is given 
by 

f(O,0) = -]og{z{6 } 0)) . (3.5) 



In particular, at /3 — 0, one obtains [18 



f(6,0) = -\og^sm( e -)) , (3.6) 



when —Ti < 9 < 7i. 
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Figure 1: U(l) Free Energy Versus 9 at (3 = 1.0 for the Naive and Binning Methods. 



We have performed Monte Carlo studies of the £7(1) gauge theory to gain insight 
into computer simulations for a system with a 9 term. The action consisted of the 
sum of the actions in Eqs. fl3.1|) and (|3.2|) . Two values of (3 were considered: (3 = 0.0 
and (3 = 1.0. Three simulation methods were employed: naive, binning and bin- 
ning with a trial probability function. Heat-bath updating was used with the naive 
and binning methods. The Metropolis algorithm was used with the trial-probability- 
binning method. The trial probability distribution Pq{Q) was chosen to be a gaussian: 
Pq{Q) oc exp(— kQ 2 ), where the constant k was appropriately selected. After ther- 
malizing the system, the number of sweeps ranged from tens of million to several 
hundred million. One sweep corresponded to updating once all the link-variables of 
the lattice. The number- crunching was done on IBM and SUN desktop workstations. 

Figure 1 plots the free energy versus 9 for (3 = 1.0, for a periodic 16 x 16 lattice. 
For the naive method, the data points correspond to short horizontal line segments. 
A total of 75 million updating sweeps were performed. The error bars were com- 
puted using a jackknife method f|6| by dividing the run into 15 data sets, each of 
which involved 5 million sweeps. The solid line is the exact analytic result. An- 
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alytic and Monte Carlo results agree for 9 less than 2.1. The agreement, which 
is excellent, cannot be seen on the scale of Figure 1. For example, at 9 = 0.5, 
f MC = 0.00242004(53) versus / exact = 0.0024196, at 9 = 1.0, f MC = 0.0096685(46) 
versus /exact = 0.0096682, at 9 = 1.5, f MC = 0.021699(61) versus / cxact = 0.021714, 
and at 9 = 2.0, juc = 0.0386(37) versus /exact = 0.0385, where the statistical un- 
certainty in the last two digits is displayed in parenthesis.^ For 9 beyond 2.1, error 
bars grew and the partition function became negative. One sees that the "barrier 
theta" 9b is about 2.1. The statistical error in Pmc(0) was 3 x 10~ 5 . Using this error 
in Eq. ( |2.10|) to estimate 9b, one finds 9b ~ 2.05. The agreement of the theoretical 
9b with the Monte Carlo value confirms the data-analysis discussion of Section II. 
One can also check Eq. ( |2.13| ) by seeing whether the ratio of the statistical error of 



P(Q) to y Pmc(Q) remains constant. Even though Pmc{Q) varied by eight orders of 
magnitude, it turned out that the ratio stayed constant to within a factor of 3. 

Figure 1 also plots the free energy versus 9 and for the binning method, again 
for j3 = 1.0 on a 16 x 16 lattice. The data points correspond to diamonds. Five 
bins of bin-size 4 were used. For each bin, 16 million sweeps were performed. Since 
5 bins were used, the statistics for this case are comparable to the statistics of the 
naive method of the previous paragraph. The 16 million sweeps were divided into 20 
sets for the jackknife analysis of errors. For 9 > 2.0, the Monte Carlo data for the 
free energy dropped below the exact result and became constant. For 9 < 1.7, the 
agreement between Monte Carlo data and the analytic result was comparable to the 
naive- method case, discussed in the previous paragraph. As one can see in Figure 1, 
the flat behavior in the free energy is anomalous for 9 > 2.1 even though the error 
bars are sizeable. At 9 = tt, the discrepancy with the exact analytic result for / is at 
the 4ct level. 

We also performed a simulation with a trial-probability-binning method, again 
with = 1.0 on a 16 x 16 lattice. The results were similar quantitatively to the 
naive case above, except that the Monte Carlo data dropped below the exact result 
at about 2.1 and then the error bars became enormous for 9 > 2.2. 

In all three of the above runs, the statistical errors in P(Q) were ordered as in 
6 For example, the 9 = 2.0 results is f M C = 0.0386 ± 0.0037. 



13 



Q 


Pmc(Q) 


-^exact 


(Q) 





1.968226(42) 


E-01 


1.968199 


E-01 


1 


1.743842(23) 


E-01 


1.743866 


E-01 


2 


1.212054(16) 


E-01 


1.212086 


E-01 


3 


6.59478(17) 


E-02 


6.59454 


E-02 


4 


2.79818(09) 


E-02 


2.79806 


E-02 


5 


9.20984(42) 


E-03 


9.20931 


E-03 


6 


2.33452(13) 


E-03 


2.33437 


E-03 



Q 


Pmc(Q) 


-^exact 


(Q) 


7 


4.51534(28) 


E-04 


4.51495 


E-04 


8 


6.58573(46) 


E-05 


6.58533 


E-05 


9 


7.13760(58) 


E-06 


7.13774 


E-06 


10 


5.64442(51) 


E-07 


5.64463 


E-07 


11 


3.18299(30) 


E-08 


3.18344 


E-08 


12 


1.24331(15) 


E-09 


1.24034 


E-09 



Table 1: P(Q) for 17(1) Model at j3 = 1.0 on a 7 X 7 Lattice. 



Eq. @. 

For the binning-method case, one can address the question of whether the flat- 
tening behavior is attributable to a value of Pmc(0) which is larger than -P(O). 
For a 16 x 16 lattice, the exact P(0) is 0.179259. The value of P MC (0) for the 
binning run was 0.179295(72). Thus, -P/v/c(0) was indeed greater than P(0). For 
the naive run, Pmc(0) = 0.179242(30), while for the trial-probability-binning run, 
Pmc(Q) — 0.179192(55), and since these two values are below the exact P(0), the 
anomalous flattening behavior is not expected to arise, in agreement with the Monte 
Carlo results. 

According to Eq. ( |2.10| ), one can increase by decreasing the size of the system. 
We therefore estimated that if the lattice size was 8, then a reliable measurement of the 
free energy could be made throughout the entire fundamental 9 region. Using a trial- 
probability-binning method on an 8 x 8 lattice, reasonable agreement between Monte 
Carlo data and analytic results was obtained in a 30-million sweep run. For example, 
at 9 = 1.0, f MC = 0.00980(5) versus / exact = 0.00967, at 9 = 2.0, f MC = 0.0382(1) 
versus / exac t = 0.0385 and at 9 = 3.0, fuc — 0.089(3) versus / exac t = 0.0814. 

Consider now the infinite- strong-coupling case (3 = 0.0. Results for this case were 
qualitatively similar to the j3 = 1.0 case: When the lattice size was sufficiently small, 
the free energy was measurable over the entire fundamental 9 region and the data 
agreed well with exact analytic calculations. When the lattice size was bigger, the free 
energy was accurately measurable only for < 9 < 9^. For example, a run, which used 
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e 




^exact 


0.2 


1.667250(80) 


E-03 


1.667220 


E-03 


0.4 


6.67569(29) 


E-03 


6.67558 


E-03 


0.6 


1.504550(71) 


E-02 


1.504530 


E-02 


0.8 


2.68107(14) 


E-02 


2.68104 


E-02 


1.0 


4.20200(25) 


E-02 


4.20195 


E-02 


1.2 


6.07374(40) 


E-02 


6.07369 


E-02 


1.4 


8.30438(62) 


E-02 


8.30436 


E-02 


1.6 


1.090390(94) 


E-01 


1.090390 


E-01 



e 


Fmc 


^exact 


1.8 


1.38842(14) 


E-01 


1.38845 


E-01 


2.0 


1.72596(23) 


E-01 


1.72603 


E-01 


2.2 


2.10467(38) 


E-01 


2.10485 


E-01 


2.4 


2.52636(72) 


E-01 


2.52674 


E-01 


2.6 


2.9917(15) 


E-01 


2.9925 


E-01 


2.8 


3.4907(35) 


E-01 


3.4926 


E-01 


3.0 


3.9436(75) 


E-01 


3.9472 


E-01 



Table 2: Measured and Exact Free Energy for (3 = 0.0 on a 4 x 4 Lattice. 



the naive method on a 30 x 30 lattice and which involved 22 million sweeps, produced 

good results only for < 9 < 0.49. The barrier value was in agreement with 9b of 

0.5, as computed from Eq. (|2.8| ) H When a binning method was used, the free-energy 

data points began to slip below the analytic results near 9 ~ 0.45. A slightly-wavy- 

but-basically-flat behavior for the free energy was observed for 0.45 < 9 < 0.80. In 

this range, the Monte Carlo data was below the exact analytic result at the 2a level. 

For 9 > 0.81, the error bars become large and the free energy was not measurable. 

When the lattice size was reduced, the free energy was measurable over a larger 9 

region. On a 7 x 7 lattice, we performed a run with the trial-probability-binning 

method. For the first 180 million sweeps, a graph qualitatively similar to the binning 

case of Figure 1 was obtained. The free energy was measurable up to 9b ~ 2.1. For 

9 > 2.1, the data points fell below the exact result and a flat free-energy curve was 

produced. From the error in P(0) of ~ 10~ 5 , the theoretically predicted value of 9b 

from Eq. (|2 . 1 ) was 2.3. When the run was continued, the flat behavior went away: 

After a billion sweeps, the data agreed with analytic results for < 9 < 2.2. Beyond 

2.2, the data fell slightly below the exact results and the error bars became large at 

9 = 2.7. For 9 > 2.7, the free energy was no longer numerically measurable. Based on 

Eq. (|2.10| ), the limiting 9b should be 2.65 for this case, in agreement with the Monte 

Carlo results. Table 1 displays the kind of accuracy with which P(Q) was measured. 

7 On this large lattice P(Q) decreased slowly for Q near zero, so that the largest statistical error 
actually occurred for P(— 4). This error was 1.4 x 1CP 5 . 
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Finally, we reduced the lattice size to 4. In a run using the binning with a probability 
distribution, the free energy was accurately obtained throughout the fundamental 9 
region. A run with 525 million sweeps was performed with 2 bins of bin-size 4. Table 
2 provides the comparison of Monte Carlo and exact results for the free energy 

The (3 = 0.0 results of this section are relevant for the CP N ^ models at /3 = 0: 
When /3 = 0, the CP^ -1 models coincide with U(l) gauge theory. A numerical 
investigation of the CP 3 model at non-zero (3 is the subject of Section IV. 

IV. The Lattice CP 3 Model 

Monte Carlo studies of the "adjoint" form of the lattice CP 3 model in the presence 
of a 9 term have been performed in ref. |28 |. For our Monte Carlo investigations, we 



have selected the "auxiliary U(l) field" formulation for two reasons: (1) it allows us 
to investigate the CP 3 model from a different-but-equally-good form of the lattice 



action, and more importantly, (2) strong coupling expansions |19| have been obtained 
for this form of the action, thereby allowing comparisons of Monte Carlo data with 
analytic results, at least for small inverse coupling (3. The CP N ~ l models without a 9 
term have been studied by computer simulations in several works. See refs. |f20[|-|[27|. 
The lattice CP N ~ l action, which we employ, is 

S = /3NJ2{4-Zx+aU{x,x + A) + c.c.) , (4.1) 

where the complex scalar fields z x satisfy Y^=\ z xi z x = 1 an d where c.c. is the complex 
conjugate of the first term in Eq. (f4.1|) . Here, the sum over A involves the d positively- 
directed nearest neighbors to x, so that A takes on the values e±, e2,. . ., e^, where 
Ci is a unit vector in the ith direction. Since we consider the two-dimensional case, 
d = 2. The field U (x, x + A) is a phase associated with link between x and x + A 
- it is the same link-variable which appears in the U(l) gauge theory of Section III. 
To Eq. ( }4.1|) , we add Sg of Eq. ( |3.2j ) to obtain the full action. Finally, we treat the 
N = 4 case, i.e., the CP 3 model. 

We used the Metropolis algorithm for all degrees of freedom. The z % x fields were 
separated into 4 real and 4 imaginary components. They were updated by performing 
rotations in each of the 28 planes of the 8-dimensional real vector space. For the U(l) 
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Figure 2: CP 3 Free Energy Versus 9 at (3 = 0.2 on 4 x 4 and 6x6 Lattices. 



fields, a binning method was used. After each time-consuming update of the z x 
fields, 10 U(l) updates, as well as Q measurements, were carried out in the fixed 
z x background. The total number of U(l) sweeps ranged from three to 150 million. 
Simulations were done for j3 = 0.2,0.6,0.7,1.0 and 1.1. The intermediate-coupling 
cross-over region is around 0.8, so that the latter two f3 values are in the weak coupling 
region where continuum scaling should set in. 

In some simulations with large volumes, the free energy was not measurable be- 
yond a certain value of 9: error bars became very large or the partition function went 
negative. The value of 9 at which this occurred was in approximate agreement with 9b 
obtained using Eq. (|2.10| ). Below, we show results for cases in which the free energy 
was measurable throughout the entire fundamental 9 region or for cases in which a 
flat behavior was observed. 

Figure 2 shows the free energy at (3 = 0.2 for 4x4 and 6x6 lattices. The 
solid line represents the tenth-order strong-coupling character expansion of ref. [ITU 



At this small value of j3, the strong-coupling expansion should be quite close to the 
exact result. On the 6x6 lattice, an anomalous flattening behavior was observed. 
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Figure 3: CP 3 Free Energy Versus 9 at (3 = 0.6 on 6 x 6 and 8x8 Lattices. 

The discrepancy between Monte Carlo data and the strong-coupling series was more 
than 10a for 9 near ir. For some reason, the jackknife error analysis produced error 
bars which did not come close to overlapping with the true results. Furthermore, 
Eq. fl2.10P predicts that 9b should be 2.1, so that measurements of the free energy 
should not be reliable for 9 > 2.1. This theoretical estimate for 9b is close to the 
point where constant free-energy behavior set in. When the lattice size was reduced 
to 4, agreement between Monte Carlo data and the strong-coupling series occurred 
throughout the entire fundamental 9 region. 

At (3 = 0.6, roughly similar behaviors in the Monte Carlo results were obtained 
except, that on the 8x8 lattice, a flat behavior in the free energy began at 9 = 1.7, 
which turned up at 9 = 2.5. The free energy was not measurable beyond 9 = 2.6. 
See Figure 3. Eq. ( |2.10| ) predicts that free energy Monte Carlo results should not be 
trusted for 9 > 1.75. For a 6 x 6 lattice, agreement with the strong-coupling series 
was obtained up to about 9 = 2.5. Beyond that point, the strong-coupling series was 
slightly below Monte data (at roughly the la level). At 9 = 2.8, the strong-coupling 
expansion of the free energy peaks and for larger 9 it decreases, while Monte Carlo 
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Figure 4: CP 3 Free Energy Versus 9 at (3 = 0.7 on a 8 x 8 Lattice. 

data continued to increase. 

At (3 = 0.7, the strong-coupling series peaks at 9 = 2.4 and then decreases signif- 
icantly. Since the free-energy should be a non- decreasing function of 9, higher-order- 
in-/3 corrections are probably important in this region, as ref. fl9| suggested might be 
the case. The Monte Carlo data began to rise above the strong-coupling results at 
around 9 = 1.2. See Figure 4. At this value of (3, which is in the intermediate coupling 
region, the Monte Carlo data is to be trusted over the strong-coupling expansion. 

At (3 = 1.0, which is on the weak coupling side of the intermediate coupling 
region, a simulation on a 20 x 20 lattice gave good results throughout the fundamental 
9 region. Figure 5 displays the free energy versus 9. The solid curve is the fit 
f{9) = 0.0039 * (1 — cos(6')). It reproduces the data within error bars. At (3 = 1.1, 
a simulation on a 40 x 40 lattice also gave good results. The free energy function 
f{9) = 0.00125 * (1 — cos(6 1 )) reproduces the (3 = 1.1 data well throughout the 
fundamental 9 region. It is interesting that a "cosine" form fits the data for (3 > 1.0. 
Such a functional form arises from a topological gas picture 0]. For f3 < 0.8, a cosine 
form did not fit the free energy data. 
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Figure 5: CP 3 Free Energy Versus 9 at (3 = 1.0 on a 20 x 20 Lattice and Its Cosine 
Fit. 



V. Remarks 

The cosine behavior of the free energy does not agree with the large N limit result, 
which gives /large n — c/N6 2 , for some iV-independent constant c. The discrepancy 
might be due to the smallness of N, which is 4 in our case. In other words, higher- 
order 1/N corrections in the large N expansion are important. We believe that this 
explanation is likely: In the leading zeroeth-order large N expansion, no 9 depen- 
dence arises in the free energy. The first 1/N correction provides the quadratic 9 
dependence. Hence, the first-order result differs significantly from the zeroeth-order 
result for large 9 and especially near 9 = n. It is thus quite possible that higher-order 
1/N corrections contribute to the large 9 region and cause the free energy to "curve 
over" like a cosine function. This would not be the first time that a large N result 
is anomalous in the CP^ -1 models. |JT], [IB, [T7| In the zeroeth-order approximation, 
the CP^ -1 models do not confine. However, the first 1/N correction does lead to 
confinement.^, |31J Nonetheless for finite N, the CP N ~ l models have a property 
called superconfinement [|T7| in which charges are so strongly bound that they can- 
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not be separated by any non-zero distance.^ Superconfinement corresponds to an 
infinite string tension. It has been conjectured that summing the 1/N series would 
strengthen the confinement and lead to superconfinement. JTTfl Probably, quantities in 
the CP N ~ l models, which differ greatly in going from the zeroeth order to the first 
order, such as the free energy for 9 away from 9 = and the string tension, are not 
reliably computable in the large N expansion. 

An interesting open question is whether there is a phase transition at 9 = n. At 



infinitely strong coupling, this is the case.[18j In the strong coupling region, finite 
volume effects round off the cusp in the free energy when sufficiently small volumes 
are used to give reliable results near 9 = it. Hence definitive conclusions cannot be 
drawn. In the weak coupling region, the cosine function fits are suggestive that the 
phase transition is absent. 

If f{9) = c'9 2 , then, as the volume V goes to infinity, 

m - exp ( (5 n 
P(o) " exp \ uv) ' 1& ' 1J 

This follows from inverting Eq. ( |2.4| ): 

P{Q) _ /! ff gexp(iflQ)exp(-V7(fl)) 

P(0) j^gexpC-V/p)) ' { - } 

As V — > oo, a saddle point expansion of Eq. (|5.2|) becomes quite accurate. If the 
minimum of f(9) is at 9 = 0, one finds 

p (Q) ( Q 2 Q 4 \ 

T(0) =exp(^-a 2 — -a 4 — -.. j , (5.3) 

where 

1 /^(O) [/(4)(Q)]2 ^(6)( 0) 1 

a2 ~VfM(0) 2V 2 [f( 2 \0)} 3 + 4F 3 [/( 2 )(0)] 5 V 3 [fW(0)} 4+ ' 

/ (4) (0) , n( 1 x ,, A , 



where /^(0) is the nth derivative of f{9) at the origin. Equation ( |5.3| ) reduces to 
Eq. Q if / (n) (0) = for n > 2. 

8 An interesting side remark is that ref. jl9| has shown that superconfinement is lost in the 
presence of a 8 term. 
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There has been some previous discussion as to whether the Q dependence in P(Q) 
is gaussian. P?|, [RJ On one hand, Eq. ( |5 . 3| ) shows that this is the case as V — > oo. 



On the other hand, if the Q dependence in P(Q) is exactly gaussian then f(9) = c'9 2 
for 9 in the fundamental 9 region, as V — > oo. Since non-quadratic behavior for 
f(9) is usually observed, the Q 4 and higher powers of Q in Eq. fl5.3| ) are important 
for determining f(9). Since the coefficient 04 of Q 4 falls off as 1/V 3 , large volume 
systems make it difficult for Monte Carlo simulations to correctly determine the 9 
dependence of f{9) for large 9. One again arrives at the conclusion of guideline (2) 
of Section II: Small volume results are to be trusted over large volume results for the 
computation of the free energy. For (3 > 0.9, the "cosine" behavior of f(9) observed 
in the Monte Carlo investigations of the CP 3 model imply a non-gaussian behavior of 
P(Q). In fact, in our simulations of both the U(l) gauge model and of the CP 3 model, 
small derivations from gaussian behavior of P(Q) were observed for all (3 values. 

VI. Discussion 

By comparing Monte Carlo simulations of the d = 2 U(l) gauge theory, we have 
verified the guidelines given in Section II: For sufficiently large volumes, a barrier 9^ 
arises, beyond which numerical results for the free energy are unreliable. The limiting 
9f, can be computed theoretically using Eq. ( |2.1U| ) with |5P(0)| estimated to be the 
statistical error in Pmc(0). When the size of the system is reduced, reliable results are 
obtained throughout the entire fundamental region for 9. Hence, when small-volume 
results differ from large-volume results, one should trust the smaller volume results. 
In some simulations, a flat free-energy behavior is observed for 9 > 9^. Comparison 
to exact analytic results demonstrates that the flat behavior is incorrect - the true 
free energy continues to rise for 9 > 9^. The flattening effect can be attributed to the 
error in Pmc(0) dominating over the errors in other Pmc{Q)- The domination of the 
error in Pmc(O) is enhanced by Monte Carlo techniques such as the binning method 
and the use of a trial probability distribution. 

The above conclusions also hold for the "auxiliary U(l) field" formulation of the 
CP 3 model: When the volume is large, anomalous flat behavior of the free energy 
is sometimes seen. At (3 = 0.2 for the 6x6 lattice run, the flat-energy behavior is 
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definitely incorrect since a comparison can be made with a reliable analytic strong- 
coupling expansion. When a smaller lattice size is used, results for the free energy are 
more accurate and flat behavior is absent. When the inverse coupling (3 is greater than 
0.6, the strong-coupling calculation of the free energy has a peak, but the Monte Carlo 
data does not. See Figures 3 and 4. At these intermediate values of the coupling, we 
believe that the peak is an artifact of truncating the series to order ten - higher-order 
contributions are probably important. 

There is a simple physical picture of why a limiting 8b arises. Current methods 
for simulating a system with 8^0 are done using the 8 = system. There should 
be a "barrier" separating the two systems. The barrier grows exponentially with the 
volume V. When V is small or when 8 is small, the barrier does not prevent the 
8 = system from sensing the physics of the 8^0 system. However, as V gets large, 
the barrier becomes more impenetrable, and for sufficiently large 8, the Monte Carlo 
simulations do not explore the phase space of the ^-system sufficiently well to give 
reliable results. 



If one applies the guidelines in Section II to the work of ref. |28], one would 
conclude the following. The simulations of the adjoint form of the lattice CP 3 model 
in refs. [2S| found the absence of a flattening behavior in the free energy for sufficiently 
small volumes. This is typified in Figure 3 of the first of refs. p8|. Guidelines (2) 
and (6) say that smaller-volume results are to be trusted over larger-volume results. 
Hence, one would conclude that the flattening behavior is anomalous. If this is true, 
it is a result of the probability method which tends to emphasize the error in Pjvrc(O). 
If the flattening behavior is anomalous, 8 C of ref. [[28| should be identified with the 
barrier theta 8b- One test of this idea is as follows. Assuming that the statistical 
errors are approximately the same for all the runs in ref. and that Eq. (|1 . 1| ) 
holds, it follows from Eq. ( 2.1C ) that V8 2 should be approximately constant for a 
fixed value of (3. At p = 2.5, for the 28 x 28, 32 x 32, 38 x 38, and 48 x 48 lattices, 
V 1/2 8 c /ir is respectively ~ 20, ~ 21, ~ 21 and ~ 22. At f3 = 2.7, for the 56 x 56, 
64 x 64, and 72 x 72 lattices, V^OJ-k is respectively ~ 35, ~ 33 and ~ 32. Hence, 
V8 2 is approximately constant. If the dominance of the error in Pmc(0) is the source 



of the free-energy flattening behavior of refs. pSI, p4J, one puzzling question arises: 
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why did all the runs display flattening behavior - one would have expected that, in 
some runs, the errors in the free energy to become large at 9 C and/or for the partition 
function to become negative. It is possible that some other subtle systematic effect 
is playing a role. The discussion here supports explanation (ii) of the Introduction, 
but one cannot definitively rule out explanation (i), namely, that a massless or light 
mode arises for 9 sufficiently large. 
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